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Abstract 

We extend the non-Markovian quantum state diffusion (QSD) equation to open quantum systems which exhibit muhi- 
channel coupHng to a harmonic oscillator reservoir. Open quantum systems which have multi-channel reservoir cou- 
pling are those in which canonical transformation of reservoir modes cannot reduce the number of reservoir operators 
appearing in the interaction Hamiltonian to one. We show that the non-Markovian QSD equation for multi-channel 
reservoir coupling can, in some cases, lead to an exact master equation which we derive. We then derive the exact 
master equation for the three-level system in a vee-type configuration which has multi-channel reservoir coupling and 
give the analytical solution. Finally, we examine the evolution of the three-level vee-type system with generalized 
Ornstein-Uhlenbeck reservoir correlations numerically. 

Keywords: Open Quantum Systems, Decoherence, Quantum State Diffusion, Master Equation, Stochastic 
Schrodinger Equation, Non-Markovian Systems 



1. Introduction 

Master equations have served as an essential tool in the dynamical analyses of open quantum systems since their 
introduction in a quantum mechanical setting by Pauli in 1928 [1]. At their core, master equations simplify the 
complexity inherent in the dynamics of a large reservoir by recording only the influence of the reservoir on the 
evolution of the system. Master equations also have the attractive feature that their solution automatically gives the 
evolution of any system observable. In contrast, while the Heisenberg picture approach may often be less complicated, 
the results are usually applicable to only one or two Heisenberg operators [2]. When the coupling to the reservoir is 
weak and when the correlation times of the reservoir are negligible compared to the evolution time of the system, the 
coupling may be said to be Markovian, and the so-called Lindblad master equation may be used to approximate the 
system evolution [3, 4]. In cases where the coupling to the reservoir is non-Markovian, alternative techniques must be 
used to determine the system dynamics. 

Methods for determining the dynamics of an open quantum system with non-Markovian coupling to a reservoir 
include the quantum trajectory approach. Quantum trajectory methods use a stochastic Schrodinger equation to evolve 
pure system states subject to a stochastic process. The evolution of the system density operator is given by the mean of 
the solution to the stochastic equation. The average of repeated numerical implementations of the stochastic evolution 
equation can be used to determine the evolution of the system when a formal solution to the stochastic equation cannot 
be found. 

A primary example of a quantum trajectories approach is the non-Markovian quantum state diffusion (QSD) equa- 
tion. First introduced by Diosi and Strunz [5], the non-Markovian quantum state diffusion equation is a generalization 
of the Markovian QSD equation [6-8]. In contrast to the Markovian QSD equation which uses a white noise process 
in the stochastic evolution equation, the non-Markovian QSD equation is built around a colored noise process. The 
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non-Markovian QSD equation has been used to examine a wide variety of non-Markovian open quantum systems, 
including the spin-N system, the N-cavity model, the N-qubit model, and many others [9-16]. The non-Markovian 
QSD equation has also been used as an analytical tool, since, in some cases, the non-Markovian QSD equation can be 
used to derive an exact master equation. For instance, the exact master equation for quantum Brownian motion [17- 
19] has also been derived using the non-Markovian QSD equation [13]. The non-Markovian quantum state diffusion 
equation is exact and applies to open quantum systems in which, a) the initial state of the total system is factorizable 
into an initial system state and a thermal reservoir state, and b) the system couples to the reservoir though a single 
channel. 

Single-channel reservoir coupling occurs when the reservoir operators entering into the interaction Hamiltonian 
are mutually proportional. Consider, for example, the following interaction Hamiltonian, 

//,„t = 2(L,„®B;';, + h.c.), (1) 

m 

where L„, are operators of the system, B,„ = YikSmkak are the reservoir operators expressed in terms of Ok, the 
annihilation operators for mode ^ of a boson reservoir, and g,,,^. is the coupling constant between system operator 
and reservoir mode k. If B,,, = Ki„B for all m, then the reservoir operators are mutually proportional and the interaction 
Hamiltonian may be reduced to interactions between a single reservoir operator, B, and an effective system operator 

Hint = ieif ® B' + h.c. (2) 



In the context of open quantum systems where the system consists of multiple particles which are coupled to a common 
reservoir, interactions of the form of (2) are often described as having collective decoherence and are an essential part 
of the theory of decoherence free subspaces [20, 21]. 

The standard non-Markovian quantum state diffusion equation may be used for interaction Hamiltonians which 
take the form of (2). Such interactions may be called single-channel reservoir coupling interactions because a canon- 
ical transformation can always be made to a new set of reservoir modes 

bj^^Ujuak, (3) 

k 

for which = V^, where rj - 1/ VSTiitF- Consequently, we have Hint - (^eff ® + h.c.) where Leff = ieff/'?- In 
contrast, multi-channel reservoir coupling occurs when a canonical transformation of reservoir modes cannot reduce 
the number of reservoir operators appearing in the interaction Hamiltonian to one (up to the Hermitian conjugate 
operation). Alternatively, we may say that the interaction Hamiltonian exhibits multi-channel reservoir coupling if the 
reservoir operators are not mutually proportional. 

Apart from the fundamental interest in solving open quantum systems with general multi-channel coupling, there is 
one system currently under both theoretical and experimental investigation which can exhibit multi-channel reservoir 
coupling. Donor-based charge quantum bits in a semiconductor host have single-channel reservoir coupling only if 
the Bohr radii of the s-wave orbitals are exactly equal for all donors [22]. Donor-based charge quantum bits are 
of particular interest due to their potential use within a scalable architecture for a quantum computer. Since multi- 
channel reservoir coupling can eliminate decoherence free subspaces which have been proposed as a strategy for 
avoiding errors in quantum computing [21, 23, 24], studying decoherence in the presence of multi-channel reservoir 
coupling is of particular importance for this system. 

In this article we will extend the quantum state diffusion method to multi-channel reservoir couplings. As is the 
case for the standard non-Markovian QSD equation, we will show that an exact master equation may sometimes be 
derived from the non-Markovian QSD equation for multi-channel reservoir coupling. We will then derive the exact 
master equation and give the analytical solution for the three-level atom in a vee-type configuration. Finally, we will 
provide a toy model for the noise correlations which will be used to examine the way in which multi-channel reservoir 
coupling can affect system evolution. 
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2. Non-Markovian QSD equation for multi-channel reservoir coupling 



We first derive the non-Markovian QSD equation associated with the interaction Hamiltonian (1). We will take an 
approach similar to the derivation of Strunz and Yu [13], and consider the case where the reservoir is initially in the 
vacuum state. A Bogoliubov transformation can be used to transform the thermal-state case to the vacuum-state case 
[12]. Adding the system and reservoir Hamiltonians to (1) we have. 



^tot — ^sys + ^int + 



Clk- 



(4) 



We go to the interaction picture of the reservoir, and write the Schrodinger equation for the total system. 



mk 



(5) 



where we use the abbreviation, = \^(t)}. We then take the component of the reservoir in the Bargmann state, 
(z|*P,), where |z> = \{zk]} = ®*z^|VAC) and Zk = e^'"'. 



" \ Jo oz„j 



In (6), we have defined z*,,, = -i Zk 8mk4^'""' 



(6) 



(7) 



and the unnormalized system state s (z| with shorthand iffiiz*) = \ilr,(z*)}. In tl/,{z*) and in the following, z* 

(z) is used as shorthand for the vector of complex Bargmann coefficients {z^j ({zk})- The functional derivative in (6) 
arises from applying the chain rule to the derivative of the Bargmann state. 



{z\ak = -^{z\ = [dsV ^^{z\- 
dz^ Jo V ^^k <s 



(8) 



As written, (6) is the exact evolution equation for \l/t{z*) and does not require interpretation as a stochastic equation, 
i.e. given an initial state of the total system I^Po), and a Bargmann state \z), (6) simply gives the evolution of the 
Bargmann component of the total system, \\l/,{z*)). Since the Bargmann state is an over-complete basis, it can be used 
to evaluate the trace over the reservoir, giving the reduced density operator of the system. 



d\p{z){z\^,){'i',\z). 



where d^z - d^z\d^Z2 ■ ■ ■ and 



(9) 



(10) 



The normalizing function p{z) appearing in (9) is of the precise form of a distribution of independent complex Gaus- 
sian random variables. We can therefore considers to be a multivariate complex Gaussian random variable, where we 
define the statistical mean M{-} of a stochastic system operator <l>(z). 



Mm 



z p(z)<l>(z). 



(11) 



In this context we interpret (z*,,) as a set of continuous Gaussian random processes with correlation functions 
0!mn{t, s) = Ai{ZmiZ*,^}, already worked out in (7). Since the noise processes all arise from the same reservoir modes. 



the cross-correlation of the noise terms given in (7) is non-zero when m ^ n (in contrast to the noise terms which 
arrive in deahng with a thermal bath [12]). Additionally, we may now interpret (6) as a stochastic evolution equation 
for \iffi{z*)), which we now call the (unnormalized) stochastic state vector Finally, the evolution of the system is found 
by taking the statistical mean of the outer product of the stochastic state vector [13], 

p, = Mmzm,(z)\]. (12) 

As in the standard non-Markovian QSD derivation we consider solutions to (6) which satisfy the following condi- 
tion: The functional derivative in the non-Markovian QSD equation may be written in terms of an operator, commonly 
called the O-operator, which depends on the Bargmann states, 

^0,n{t,S,Z*)>^,. (13) 

where 0„,(s, s,z*) = L„i so that (6) is able to reproduce the QSD equation in the Markov limit [6-8]. 
By defining 



QJt, z') = j^ds ^ a,„„(t, s)0„{U s, z*)j , (14) 

we can write the non-Markovian QSD equation under the ansatz, 

ihdiff,^H,s(t,z*)ifr„ (15) 
where Hgff{t, z*) is the effective stochastic Hamiltonian, 

//eff(f, z*) = //sys + ih ^ (z,;,L„, - LlQJt, z*)) . (16) 

nt 

To be self-consistent, the solution condition specified in (13) must satisfy the following conditions, 

d,-^ = T^^'-A" (17) 
giving rise to the consistency equations for the O-operator 0,„{t, s, z*), 

d,0„, = [H,dt,z*), 0„,] - Yj i^l4rQnit,z). (18) 

Finally, by making a Girsanov transformation [10, 25], we can write the nonlinear non-Markovian QSD equation, 
which is more suitable for numerical simulation, 

d,i(r, = --^//sysfAr + [2"„A(i„,) - A (A(L;'„)e„,(f, ?)) ]>J^„ (19) 

m 

where A(A) - A - (A) for any operator A and z*„, = z*„ + 2,, amnit, s){L'i„)ds. 

Equations (15-16) and (18-19) are the main results of this article. They allow systems with multi-channel non- 
Markovian reservoir coupling to be evaluated numerically using finite-difference methods. Additionally, in some 
cases, the quantum state diffusion equations allow for more efficient simulation as compared to master equations in 
cases where the master equation may also be derived [9]. 

When the O-operators are noise-independent, a master equation may be derived, as is also the case for single- 
channel coupling. Following the treatment of [13] we evaluate the time-derivative of the density operator in (12) with 
the help of (15) to arrive at 

|a = -{[H.y.,p,] + J] {]^M{Q„,{t,z')P,],Ll] + h.c), (20) 
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Figure 1: a) Thi'ee-level system in vee-type configuration, b) Multi-level vee-lambda system with 3 upper levels and 4 lower levels. The lines 
connecting levels indicate possible transitions. 



where f , — \il/t(z*)}{tf't(z*)\ and we have used the Novikov theorem and the operator ansatz (13) [13, 26] to write 



When the O-operators are noise independent, then Q„{t, z") - QmiO and we can write an exact master equation. 



which is the multi-channel analog of the exact master equation found by Strunz and Yu [13] for single-channel reser- 
voir couplings. 

3. Three-level vee-type system 

We now apply the previous results to demonstrate the utility of the QSD approach for open quantum systems with 
multi-channel reservoir coupling. We consider the three-level system in a vee-type configuration with multi-channel 
reservoir coupling, as shown in Fig. la). We consider the case where the decoherence is generated by the lowering 
operators, L,„ = \3){m\. 



where the excited levels |1) and |2) are separated from the ground state |3) by energies fi-wi and ha>2, respectively, so 
that i/sys - |l)(l| + fi.W2|2)(2|. We note that in the following, summations will always be over indices corresponding 
to the two upper levels, 1 and 2. When we expand the O-operator as 



(21) 




(22) 



m 



Hint = 13X11 ®B|" + |3)<2|®BT+h.c., 



(23) 




(24) 



p 



we find that the consistency equations reduce to 




(25) 



m 



where fqp(s, s) = dqp, and we have defined 




n 



(26) 



so that Q„,{t, z') - Yjp Fmp{t)Lp. The resulting master equation is given by 



(27) 



mp 
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With only slight adjustment, the O-operator in (24) can be generalized successfully to any system which can be 
divided into upper and lower levels with transitions forbidden within the upper and lower levels, respectively, and in 
which the decoherence generators are lowering operators between upper and lower levels [27]. We call such systems 
multi-level vee-lambda systems since they are a generalization of both the vee and lambda systems. Figure lb gives 
an example of a 7-level vee-lambda system. The resulting master equation will be of the exact form of (27), albeit 
with the summation running over all upper and lower levels. 

Though the master equation in (27) appears similar in nature to approximate master equations (such as the Born- 
Markov or Lindblad master equations), (27) is exact. The time-dependence present in the coefficients encapsulates 
the total effect of the reservoir, and, in particular, any non-Markovian effects. 

For the three-level vee-type system we can find an analytic solution to the master equation. We first define the 
operator 

F(t)^Yj^r„„{t)\m){nl (28) 

m 

which allows us to write 



5,p = - [H,y,,p] - (Fp + pF^ ) + tr I + pF' \ |3)<3|. (29) 



The solution to (29) with initial condition p(0) = po is given by 

pit) - p(t)pM + (l - p(f))|3><3| (30a) 

where 

pit)^tT[olO,po\, (30b) 
OipoO] 

Peit) = r V r (30c) 

tT[OlOrPo\ 



and 



Of - exp 



-jH,,,t- \ dsF{s) 



- fd 

Jo 



(30d) 



Equation (30) gives the general solution of a vee-type three-level system coupled to a vacuum harmonic os- 
cillator reservoir under the rotating wave approximation. General properties of vee-type systems can be deduced 
from the form of (30). Since H^ys and F{t) have support in only the system subspace defined by the excited states, 
"He = Span[{|l), |2))], if the system initial state is excited, po e , thenp(.(f) is also an excited state, P(,(f) € TYp. Addi- 
tionally, if the initial state is pure, po = |iAo)(iAoL then pp(f) = \4it){4't\ is also pure, where oc (9,|i^o)- Consequently, 
we can interpret the solution as follows; at time t an initially pure and initially excited state will have either decayed 
to the ground state with probability (1 - pit)) or, with probability p{t), have undergone unitary evolution within the 
excited state subspace. 

4. Ornstein-Uhlenbeck Correlations 

We now consider a toy model for the noise correlations which will allow us to examine the effect that non- 
Markovian and/or multi-channel reservoir coupling has on the system dynamics. The noise model we will investigate 
is the generalized Ornstein-Uhlenbeck [OU] correlation function, 

a,™(r) = (.-<'''"-"'"'^0(T) + .<''"-'"">^0(-T)) , (31) 
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where t — t - s, and 6{t) is the Heaviside function. The generahzed OU correlation function reproduces the standard 
OU correlation function when n - m, 

«„™(t) = |^„,|2^e-'''"We-'""^ (32) 

We define Y„, - \k„,'^ since it corresponds to the decay rate in the standard OU correlation function. D.„, is the central 
frequency of the OU correlation function. We also see that 1 Ijm is the correlation time of Bm, as y,„ — > oo the standard 
OU correlation function goes to a delta function ai„„,(T) — » r„,S(T). For this reason a finite value of y,„ signals the 
departure from the Markov limit of the reservoir coupling to level m. 

The parameters in (31) can also be described in terms of a quasi-Lorentzian model for the coupling coefficient, 

gm((^) = . ^ ^ —, (33) 

where «'(w) is the spectral density of reservoir modes, /<•„, is a coupling strength (with units of VHz), and y,„,Qm 
(where jm, Qm > and have units of Hz) are the bandwidth and central frequency of the reservoir coupling. By letting 
the summation over k in (7) go to an integral over all frequencies (positive and negative), the Lorentzian coupling 
coefficient (33) reproduces the noise correlations in (31). In the following we will describe the parameters of the OU 
correlations in terms of their properties in the coupling coefficients. 

Using the generalized OU correlation function, the evolution of the time dependent operator F{s) in (30) is given 

by 

(34) 

/' 

This set of coupled non-linear first order equations can be solved numerically for particular choices of the six param- 
eters in the correlation function {*■„,, y„,, Q,„) and the two atomic energies {w,,,). The solution can then be combined 
with the analytical solution in (30) to give the complete dynamics of the atomic density operator 



5. Single-channel reservoir coupling 

We now investigate some choices for the coupling parameters to show how multi-channel reservoir coupling 
manifests itself in this system. To do this we first examine the reduction of the consistency equation to the single- 
channel coupling case, B,„ - k,„B, which implies that y,„ = y and Q,„ = Q. If, in addition, we consider the case where 
levels 1 and 2 are degenerate, lo,„ = wo, then we can derive the evolution of the coefficients F„,„{t) analytically. By 
defining A - (wq - fi),yS = -(y - iA)/2, and jf = K~y/2 -ff (where k~ - \ki\^ + \k2\^ = Ti + Fa), it can be shown that 

F,„n(t) = ^2(0, (35) 
where Q(t) satisfies the differential equation 

d,Q^(Q+/3)- + T^-, (36) 
with solution 

y / sin 77f \ 

Q(t)^K^7:\ —. , (37) 

2 \ 7/ COS rjt - f) sm jjt ) 

under the initial condition 2(0) = 0. 

Q{t) governs the decoherence behavior for many qualitatively different parameter regimes. For some of these 
regimes, Q{t) exhibits complex infinities similar to those described by Diosi et al. [9]. These infinities pose no 
difficulty to the evolution of the quantum system as discussed by Strunz etal. [11]. That conclusion is confirmed here 
since 



exp 



Jo 



dsQis) 



— (ri cos r]t - /3 sin T]t), (38) 
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Figure 2: (a-d) Plot of P33 (solid), pn (dashed), p22 (dotted), and |pi2l (dot-dashed) when r„, = wq, A = tJo/lOO, and y,,, = y = {0.1, 0.2, 1, 5)(Jo 
for the initial state p(0) = |1)(1|. In all plots the system state in the long-time limit is given by pn = P22 = IP12I = 0.25 and P33 = 0.5. 



SO that infinities in Q(t) correspond to zeros in the evolution of p,. 
If we define the state 



K ^ 

m 

we find that we can write the master equation as 

d,p, = -UH,y„p,] + Q{t)[[L^ffp„Llff] + [Leff,Aieff]), 



(39) 



(40) 



where Lgif - |3){0^|. Written this way it is easy to see that \<p^} cc k2\1}-ki\2} is a trivial solution to (40). Alternatively, 
it can be shown that satisfies the conditions of a decoherence free subspace [20, 23, 24]. Additionally, since any 
excited state will eventually decay to a mixture of and the ground state, the vacuum reservoir induces coherence 
between levels 1 and 2, for all values of y > 0. This effect was first discovered for the three-level vee-type system 
by Agarwal [28] in the Markov limit and is known as vacuum-induced coherence'. It should be noted that the three- 
level vee-type system under a rotating wave approximation is essentially equivalent to the problem of 2 two-level 
atoms sharing a single excitation under a rotating wave and an essential states approximation. Thus, vacuum-induced 
coherence in the degenerate three-level vee-type system is directly analogous to the super- and sub-radiant effects in 
degenerate two-level atoms first emphasized by Dicke [30]. 
In the Markov limit, when 7 » coq, Q, k^, we find that 



Qit) - 2 ' 



so that the master equation in the Markov limit is given by 
dtPt 



My.,p,] + ^^{[L eSPl, L^g] + [Le:S,p,L^^gfj , 



(41) 



(42) 



as can be verified using traditional methods [31]. We see in (42) that the decay rate out of is the sum of the 
individual decay rates out of levels 1 and 2. 

We now consider the case where the decay rates are equal T,,, = T. Under this condition, the density matrix 
equations in the Markov limit (42) are exactly those found by Agarwal [28]. We plot the evolution of the level 
populations, pn, P22, P33, and the excited state coherence |pi2| in Fig. 2a-d for various choices of the width of the 
coupling coefficients, 7, for the initial state p(0) = 

We choose the coupling to be very strong, F = a>o, so that the decoherence evolution takes place in a few hundred 
cycles as opposed to many thousand cycles when F «: ojq. We also choose the central frequencies of the coupling 
coefficients to be nearly resonant with the excited state energy, A = (cjo-Q.) - -wo/lOO. As the width of the coupling 



' In the literature there appear to be two scenarios related to vacuum-induced coherence. The other one includes classical fields which participate 
in the effect. See, for instance Ref. [29]. 
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Figure 3: (a-d) Plot of P33 (solid), pu (dashed), P22 (dotted), and |pi2| (dot-dashed) when F,,, = wo, A = (Jq/IOO, yi = 5(t)o, and 
y2 = |5,2.5, l,0.1|wo for the initial state p(0) = |</i")(0"|. In (a), P22 = Pii = IP12I = 0.5 as indicated by the overlap of the dashed, dotted, 
and dot-dashed lines. 

coefficient increases from y - 0.1 wo we see that the highly non-Markovian evolution in Fig. 2a is replaced by the 
Markovian decay in Fig. 2d where y - Scoq. We also see that after a long time the total system state evolves to a 
mixture of\(p^) and the ground state, 

P, + \\^){% (43) 

but that the time required to arrive at the long-time limit depends upon the correlation time, 1 /y, of the reservoir 
operator B. 

6. Multi-channel coupling with variable widths 

We now consider the case where the widths of the correlation functions begin to deviate from one another, but all 
other parameters are held fixed. When the widths of the coupling coefficients are not equal, the reservoir coupling is 
no longer to a single channel. We investigate the case where the width of the coupling to level 1 is given by yi - Soiq 
and the initial state of the atom is |0 ) to examine how multi-channel coupling affects the evolution of the decoherence 
free state. In Fig. 3a-d we plot the level populations p,„„,, and the coherence between excited states |pi2| as the width 
of the second coefficient, 72, decreases from the width of the first coefficient. As already discussed, the initial state 
satisfies the conditions of a decoherence free subspace when yi — y2, as shown in Fig. 3a. As the width of the second 
coefficient decreases from the first, the system begins to decay to the ground state. When the width of the second 
coupling coefficient becomes very narrow relative to the first, 72 = O.lyi, the decay of pn roughly reduces to what it 
would be if the decay from level 2 was forbidden, r2 = 0. Consequently, we see that a common noise source mitigates 
the decay from level 1, even when a decoherence free subspace is not supported by the interaction. 

7. Multi-channel coupUng with variable central frequencies 

Finally, we investigate the case where the width of the coupling coefficients are equal y„, = y, but the central 
frequencies are allowed to shift. When the central frequencies are not identical, the reservoir coupling is no longer to 
a single channel. We examine the case where the central frequency of the second coupling coefficient, D.2, decreases 
from the first which remains near resonance, Qi - l.Olwo, and plot the results in Fig. 4a-f. Since the shift of the 
coupling coefficient relative to the width will be important, we investigate two cases: In Fig. 4a-c we let y - and 
let Q.2 - {.67, .33, 0)^0, and in Fig. 4d-f, we halve the bandwidth, y - O.Swo and let the central frequency have the 
same values as before, Q.2 = {.67, .33, 0)wo- 

Regardless of the width of the reservoir coupling, when central frequencies are the same, the initial state |0 ) does 
not decay, as shown for y = Swq in Fig. 3a. As the central frequency of the second channel deviates from the first, the 
frequency separation relative to width becomes important. In Fig. 4a the frequency separation IQ2 - Qi | is only a third 
of the width. In contrast. Fig. 4d has a frequency separation which is 2/3 of the width. Consequently, the erosion of 
the decoherence free state in Fig. 4d is more severe than in Fig. 4a. This property is also apparent in comparing the 
rate of decay in Fig. 4b and e, and c and f. 
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Figure 4: (a-f) Plot of P33 (solid), pn (dashed), pii (dotted), and \p\i\ (dot-dashed) when T = Ci>(,, and Cli = l.Olwo for the initial state p(0) = 
\(p~)(ip~\. Plots (a-c) in the left column have y = a)o whereas plots (d-f) in the right column have y = 0.5(t<o- The plots in the top (a,d), middle (b,e), 
and bottom (c,f) rows have Q.2 = {.67, .33, 0)(i;o, respectively. In all plots, pi 1, P22. and Ipnl decay in nearly the same fashion from 0.5 as indicated 
by the overlap of the dashed, dotted, and dot-dashed lines. 



When the separation of the central frequencies is greater than their width, increasing the separation does not have 
a strong additional effect on the system evolution. This can be seen in Fig. 4d-f; the rate of decay in Fig. 4e and 
f remains roughly the same even though the frequency separation has increased in Fig. 4f by the same amount as 
between Fig. 4d and e. 

Figs. 2-4 are also consistent with the result that the decoherence free subspaces are robust. Bacon et al. [32] 
showed that perturbations to interactions which support decoherence free subspaces do not degrade the decoherence 
free subspace to first order in the perturbation and to first order in time. To see how this applies to the plots above, we 
rewrite the interaction Hamiltonian as 

//.nt^Z/rm+^rnt (44) 

where Hf^^ - (\3}{(p'^\ ig) b\_+ h.c), and where B± cc Bi± 82- When Bi ^ B2, then B ^ 0, and we can consider H^^^ as 
a perturbation to H^^^. Consequently, to the degree that Bi ^ B2 we should expect similarly reduced decay out of 
which is a decoherence free subspace of H^^^. This is verified in Fig. 3b, and Fig. 4a, b, and d. 



8. Conclusion 

We have extended the non-Markovian quantum state diffusion equation to open quantum systems with multi- 
channel reservoir coupling. Multi-channel reservoir coupling occurs when a canonical transformation of reservoir 
modes cannot reduce the number of reservoir operators appearing in the interaction Hamiltonian to one. Apart from 
fundamental interest in solving the most general type of reservoir coupling, open quantum systems which can exhibit 
multi-channel reservoir coupling are now under both theoretical and experimental investigation [22]. Additionally, 
understanding the effect that multi-channel coupling has on decoherence will likely be important for developing 
practical quantum information systems. 

For open quantum systems with multi-channel reservoir coupling which admit a noise-free O-operator, we have 
derived the exact master equation. We then considered non-Markovian evolution in the three-level vee-type system, 
finding the exact master equation in terms of time-dependent coefficients. By reformulating the master equation in 
terms of a time-dependent decay operator, we found the analytical solution to the general master equation for the 
three-level vee-type system with multi-channel coupling to a harmonic oscillator reservoir 

Using generalized Ornstein-Uhlenbeck noise we demonstrated how the solution to the master equation for the 
three-level vee-type system may be used once the noise correlation function has been determined. When the correla- 
tions are identical (up to a constant), the multi-channel reservoir coupling reduces to single-channel reservoir coupling. 
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and the coefficient equations (34) can been solved analytically. In this limit, the decoherence free state correspond- 
ing to vacuum-induced coherence is supported by the interaction. As the spectral widths of the Ornstein-Uhlenbeck 
correlations deviate from each other, or as their central frequencies separate, the decoherence free state decays. When 
the correlations deviate strongly, in width or in central frequency, the decay proceeds much more rapidly than when 
the correlations are nearly the same. This demonstrates the robustness of decoherence free subspaces for Ornstein- 
Uhlenbeck noise correlations. 
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